arXiv:1502.03099v2 [astro-ph.CO] 23 Sep 2015 


YITP-15-7 


Cosmology with anisotropic galaxy clustering from the combination of power 

spectrum and bispectrum 

Yong-Seon Song^, Atsushi Taruya^’^, Akira Oka^ 

^Korea Astronomy and Space Science Institute, 

Daejeon 305-348, R. Korea and 
^ Yukawa Institute for Theoretical Physics, 

Kyoto University, Kyoto 606-8502, Japan 
^ Kavli Institute for the Physics and Mathematics of the Universe, 

Todai Institutes for Advanced Study, 
the University of Tokyo, 

Kashiwa, Chiba 277-8583, Japan 
* Department of Physics, 

University of Tokyo, Tokyo 113-0033, Japan 

(Dated: September 25, 2015) 

The apparent anisotropies of the galaxy clustering in observable redshift space provide a unique 
opportunity to simultaneously probe cosmic expansion and gravity on cosmological scales via the 
Alcock-Paczynski effect and redshift-space distortions. While the improved theoretical models have 
been proposed and developed to describe the apparent anisotropic clustering at weakly non-linear 
scales, the applicability of these models is still limited in the presence of the non-perturbative 
smearing effect caused by the randomness of the relative velocities. Although the cosmological 
constraint from the anisotropic clustering will be improved with a more elaborate theoretical model, 
here we consider an alternative approach using the statistical power of both the power spectrum and 
bispectrum at large scales. Based on the Fisher matrix analysis, we estimate the benefit of combining 
the power spectrum and bispectrum, finding that for the future spectroscopy survey DESI (Dark 
Energy Spectroscopy Instrument), the constraints on the cosmic expansion and growth of structure 
will be improved by a factor of two. This approach compensates for the loss of constraining power, 
using the power spectrum alone, due to the randomness of the relative velocities. 

PACS numbers: 


I. INTRODUCTION 

In our current understanding of the universe, an un¬ 
known substance called dark matter dominates over the 
standard model particles at the present epoch. Despite 
many theoretical and observational efforts, the origin of 
dark matter is not yet clarified. Also, the existence of 
dark energy, which is supposed to drive the cosmic ac¬ 
celeration, indicates our incomplete understanding of the 
gravity on cosmological scales 0 . It may imply mod¬ 
ifications to Einstein’s theory of General Relativity. A 
further insight into the origin and nature of dark energy 
or validity of general relativity is essential, and this is 
one of the primary goals in next-generation cosmology. 

The large-scale structure offers an opportunity to 
probe these issues by looking at the anisotropic galaxy 
clustering in redshift space The observed galaxy 

distribution via the spectroscopic measurements is ap¬ 
parently distorted due to the peculiar velocity of galax¬ 
ies along the line-of-sight direction, referred to as the 
redshift-space distortions (RSD). While the RSD compli¬ 
cates the interpretation of the small-scale galaxy cluster¬ 
ing, on large scales, the strength of anisotropies is simply 
characterized by the linear growth rate / = dIn/dina 
P, providing us a unique opportunity to probe the 
growth of structure, where dm and a are the linear den¬ 
sity field and scale factor of the Universe, respectively. 


On the other hand, the large-scale galaxy clustering data 
imprints a fossil record of the primeval baryon-photon 
fluid around the last-scattering surface, known as the 
baryon acoustic oscillations (BAO) [§-[ll|. The charac¬ 
teristic scale of the BAO can be used as a standard ruler, 
which enables us to determine the geometric distances of 
high-z galaxies with a greater precision. The key point to 
determine the geometric distances is to measure the clus¬ 
tering anisotropies over the BAO scales. Notice the fact 
that the anisotropies of the clustering pattern also arises 
from the apparent mismatch of the underlying cosmologi¬ 
cal model when we convert the redshift and angular posi¬ 
tion of each galaxy to the co-moving radial and transverse 
distances. This is the so-called Alcock-Paczynski (A-P) 
effect [l2| , and with a prior knowledge of the characteris¬ 
tic scale of the BAO, the Hubble parameter H{z) and an¬ 
gular diameter distance Da{z) of the high-z galaxies can 
be separately measured. Thus, the anisotropic galaxy 
clustering can serve as a dual cosmological probe from 
which we can explore the origin of cosmic acceleration 
from the viewpoint of both dark energy and modification 
of gravity. 

In order to simultaneously extract information on both 
the growth of structure and cosmic expansion, a detailed 
theoretical model of the anisotropic clustering is crucial. 
While the scales of our interest are rather close to the 
linear regime of structure formation, there appear small 
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but non-negligible nonlinear systematics that must be 
corrected or subtracted, and these including the nonlin¬ 
ear gravitational clustering. This is one of the main rea¬ 
sons why perturbation theory has become very popular 
recently as a theoretical template of the power spectrum 
or correlation function beyond the linear theory (e.g., 
[l3l - IHj| l. Still, however, perturbation theory has its own 
limitations, and we cannot apply it to the small scales 
beyond the weakly nonlinear regime. Furthermore, the 
RSD described by the nonlinear mapping from the real 
to redshift spaces, gives an additional complication that 
leads to the non-trivial cross talk between small- and 
large-scale clustering. As a result, even at large scales, 
the clustering amplitude is significantly reduced along the 
line of sight, known as the Finger-of-God effect. While 
several treatments, based on perturbation theory, have 
been proposed in order to precisely model the nonlin¬ 
ear RSD (e.g., [12-113), the Finger-of-God effect, being 
most likely ascribed to the virialized random motion the 
galaxies inside halos, prevents us from a detailed mod¬ 
eling within perturbation theory. Hence, a phenomeno¬ 
logical description characterizing the Finger-of-God sup¬ 
pression needs to be introduced, and in order to avoid 
any unwanted systematics, we have to conservatively re¬ 
strict the cosmological data analysis to the larger scales, 
k < 0.1/iMpc“^, for instance. This, needless to say, sig¬ 
nificantly reduces the statistical power to constrain dark 
energy or to test general relativity. 

In the persuit of extracting maximal cosmological in¬ 
formation from the anisotropic galaxy clustering data, a 
simple but potentially powerful approach is to make use 
of the benefit of combining both the power spectrum and 
bispectrum on large scales. Although the initial condi¬ 
tion for perturbations is supposed to be nearly Gaussian, 
non-vanishing bisepctrum naturally arises from the non¬ 
linear mode coupling through the late-time gravitational 
evolution. At the weakly nonlinear regime, the bispec¬ 
trum still contains statistical information similar to the 
power spectrum, and thus the BAG feature should be 
clearly manifest. Hence, the bispectrum in redshift space, 
as an actual observable, can be used as an alternative 
tool to constrain the geometric distances and growth of 
structure via the A-P effect and RSD. Thus, in combin¬ 
ing the bispectrum and the power spectrum, we expect 
a substantial on the resultant cosmological constraints. 
Indeed, in the context of cosmology with galaxy redshift 
surveys, the benefit of using the bispectrum in cosmologi¬ 
cal data analysis has already been studied in several work 
(e.g., [sil - l^ l. Among these, Sefusatti, Crocce, Pueblas 
& Scoccimarro [ 11 ] focused on the angle-averaged bispec¬ 
trum in redshift space, and discussed a potential power 
of bispectrum to constrain multiple cosmological parame¬ 
ters, ignoring A-P effect. On the other hand, taking fully 
account of both the A-P effect and RSD, Greig, Komatsu 
& Whyithe [Hj considered cosmology with Lya emitting 
galaxies, and specifically studied the impact of radiative 
transfer effects on the observed clustering of Lya emit¬ 
ting galaxies. It has been shown that the bispectrum 


is helpful to distinguish between gravitational and non- 
gravitational effects, thus breaking the parameter degen¬ 
eracies. 

In this paper, we present the combined results of the 
power spectrum and bispectrum to constrain the geomet¬ 
ric distances and growth of structure. The present paper 
is partly similar to [slj - but is rather different in sev¬ 
eral aspects. To be specific, we consider a future accessi¬ 
ble spectroscopic survey like Dark Energy Spectroscopy 
Instrument (hereafter DESI), which would be the best 
suited to probe the cosmic acceleration around z ^ 1. 
Gonsidering this survey setup, we discuss the impact of 
the Finger-of-God effect on the estimation of cosmolog- 
ical parameters, which has not been considered in [33|| . 
We show that while the uncertainty of the Finger-of- 
God effect in power spectrum is mostly degenerate with 
the coherent motion as a probe of the growth of struc¬ 
ture, the combination of power spectrum and bispectrum 
breaks this degeneracy, thus improving the measurement 
accuracy of the coherent motion by a factor of two. As 
for the constraints on geometric distances, substantial 
improvement is found, consistent with previous works. 
Furthermore, the role of the non-vanishing cross covari¬ 
ance between power spectrum and bispectrum is studied, 
and the relative impact in estimating statistical errors is 
quantified, finding that the influence of cross covariance 
is small enough in the weakly nonlinear regime. 

This paper is organized as follows. In Sec. |IT1 the ba¬ 
sic setup of the forecast analysis is presented. Based 
on perturbation theory, theoretical models for power 
spectrum and bispectrum are given, taking account of 
both the A-P effect and RSD. Then, the basic formal¬ 
ism for Fisher matrix analysis is described, including 
the cross covariance between power spectrum and bis¬ 
pectrum. Sec. m presents the main results for Fisher 
matrix analysis. After comparing the signal-to-noise ra¬ 
tios of power spectrum with those of bispectrum, we show 
the expected constraints on geometric distances {Da and 
H) and growth of structure / in DESI-like experiments. 
The impact of Finger-of-God effect is discussed in detail, 
and role of the bispectrum is clarified. Finally, Sec. m is 
devoted to the summary and discussion. 


II. CONSTRAINTS ON COSMOLOGY USING 
POWER SPECTRUM AND BISPECTRUM 

A. Power spectrum 

On large scales of our interest, the density and ve¬ 
locity fields are basically the small perturbations to the 
homogenous background. When the higher-order con¬ 
tributions are ignorable, the power spectrum is simply 
described by the linear theory, and in redshift space, it 
is given by 

P'”(fc) = Z2(fc)P(fc), 


( 1 ) 
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where Zi{k) is defined by, 


B. Bispectrum 


Z^{k)=b + f^l\ (2) 

The linear galaxy bias, denoted by 6 , represents the en¬ 
hancement of the clustering amplitude relative to the 
mass density field (5ni. The function / is defined by 
the logarithmic derivative of the linear density field 
with respect to the logarithm of scale factor, i.e., / = 
(iln(5m/dlna. The directional vector k is decomposed 
into (fc, /i) where fi denotes the cosine of angle to the line 
of sight. 

In practice, the applicability of the linear theory ex¬ 
pression in Eq. m is restricted to a narrow range of 
scales. This is because the mapping of statistical quan¬ 
tities from real to redshift space is intrinsically nonlin¬ 
ear. Even at large scales, higher-order contributions to 
the mapping formula are not negligible. Eor this rea¬ 
son, there have been several improved models of RSD 
that have been proposed that add correction terms to 
the Kaiser formula in Eq. ©• Taking account of those 
effects is thus crucial and essential for an unbiased pa¬ 
rameter estimation in the practical data analysis. On 
the other hand, the estimation of statistical errors is not 
usually much dependent on those elaborate factorized 
formulations, because the statistical error of each param¬ 
eter mainly comes from the measurement uncertainties, 
including the cosmic variance and shot noise. Unless a 
significant contribution of higher-order corrections arises, 
the structure of parameter degeneracies will remain un¬ 
changed. Hence, in this paper, we do not consider such 
higher-order corrections. 

Nevertheless, the suppression of clustering ampli¬ 
tude due to the random motion is known as a non- 
perturbative effect, which significantly affects the power 
spectrum even at large scales, and should be accounted 
in the basic formulation of Eq. ©■ Here, we assume that 
this EoG is^vm as a factorized form, and multiplied as 


While the initial condition for perturbations is as¬ 
sumed to be Gaussian, gravitational evolution naturally 
induces mode-mode coupling, giving rise to the non¬ 
vanishing bispectrum. Eurthermore, coupled with galaxy 
bias and RSD, the bispectrum in redshift space becomes 
rather complicated. The resultant leading-order expres¬ 
sion for the bispectrum (e.g., [3l|), valid at large scales, 
is given by. 


H^^(fcl,fc2,fc3) = 2 Z2{ki,k2)Zi{ki)Zi{k2)P{ki)P{k2). 


cyclic 


( 5 ) 


The kernel Z 2 is defined as, 

Z2{ki,kj) = ^ +bF2{k^,kj) + f^j.fjG2iki,kj) ( 6 ) 

+/Mj) + ^(^ +) J 

iCj Kj 

where we define fii = (ki ■ z)/ki, kij = ki + kj, fiij = 
{kij ■ z)/kij, with z being the line-of-sight unit vec¬ 
tor. Here, we incorporate the uncertainty of the non¬ 
linear galaxy bias characterized by 62 into the kernel 
Z 2 , adopting the local bias prescription (e.g., [HI), i.e., 

Sg = b5ra + {b2/2) -. In the above, the functions F 2 

and G 2 are the standard PT kernel in real space, given 

by, 


f kij 


F2{h,kj) 

G2{h,kj) 


2 
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kj 


kj 


kj 
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( 7 ) 

( 8 ) 


with rjij = {ki ■ kj)/{kikj). Note that the configuration 
of bispectrum satisfies the triangular condition, which is 
expressed by the directional vector constraint. 


P{k) = ■ 


( 3 ) 


The Dp^Q{k) is given by the Gaussian form as, 


-Dfog(^) = exp - {kjiapf , 


( 4 ) 


where Op denotes the dispersion of the one-point PDE 
of the velocity in one-dimension. Note that at smaller 
scales, the virial motion of galaxies inside a cluster of 
galaxies also leads to a suppression of the power spec¬ 
trum in redshift space. When (kfiap)^ <C I, the leading 
order term of Eq. ([3]) is dominant over all other higher 
orders, and the estimated errors are immune from the 
exact functional form of Eq. The linear Cp is used 
for the fiducial value. 


^1 + ^2 + ^3 — 0 . (9) 

In contrast to the redshift-space power spectrum, the 
influence of nonlinear RSD on Eq. ([5]) is not yet fully 
understood and studied in detail. Although it deserves 
further investigation, we can make the best guess on the 
possible damping effect due to the random motion of 
galaxy. The EoG effect in the bispectrum is assumed 
to be Gaussian as [4lj |. 

77foG(^i42,4) = exp [-{kfjil + kjul + klnl)al] .(10) 

Then the observed bispectrum is given by, 

B{ki,k2,k3) = Dp^Q{ki, ^ 2 ,*2,4) • (H) 

Again, when {k^ap)'^ I, the detailed functional form 
of Dp^Q is not important for our estimation. 
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C. Alcock—Paczynski test 


where a(/ii) is defined by, 


In addition to the anisotropies induced by the RSD, 
the observed galaxy clustering also exhibits anisotropies 
through the Alcock-Paczynski (A-P) effect. This can 
happen if the background expansion of the real universe 
differs from the fiducial cosmology used to convert the 
redshift and angular position of each galaxy to the co¬ 
moving radial and transverse distances. 

While this effect leads to the modulation in the shape 
and amplitude of the power spectrum and bispectrum, if 
the shape of these quantities is a priori known, it offers 
a unique opportunity to measure the angular diameter 
distance Da{z) and Hubble parameter H{z) of distant 
galaxies at redshift 2 ; using the characteristic shape of 
the galaxy clustering in both the radial and transverse 
directions. Furthermore, providing information on the 
evolution of density and velocity fields, the two types of 
apparent anisotropies (i.e., RSD and A-P effects) become 
distinguishable, and the geometric distances Da and H 
can be separately and accurately determined. This is in¬ 
deed possible if we know at least the broadband shape 
of spectrum. In other words, given a accurate theoreti¬ 
cal template which describes the broadband shape of the 
power spectrum and bispectrum, the simultaneous con¬ 
straints on the geometric distances and growth of struc¬ 
ture are made possible. We dub this method as a broad¬ 
band A-P test [ 3 ^ . 

The anisotropies in the power spectrum caused by the 
A-P effect are modeled as follows. Denoting the true 
power spectrum by P , the observed power spectrum be¬ 
comes 




^true \ 

i/fid ) 


jj&d 

7~)true 


2 

^( 9 , v) , 


( 12 ) 


a(/ii) = 


' Df 




jjfid 


jjdd 

7~)true 


1/2 


The cosine of angle between two vectors, Vij = 
is given by, 


nfid 

^A 

7~)true 

^A 

jjtr. 


Vij 


fffid 




nfld 

^A 


7~)true 

^A 






■ (16) 


Here, we define rjij = [ki ■ kj)/(kikj). 


D. Covariance matrix and Fisher matrix analysis 

In this paper, to elucidate the potential power of the 
bispectrum to constrain cosmology, we shall specifically 
consider DESI as a representative future galaxy survey, 
and proceed to the Fisher matrix analysis. The primary 
science goal of the DESI experiment is to clarify the na¬ 
ture of dark energy and/or gravity through the A-P and 
RSD effects, and starting in 2018, it will obtain optical 
spectra for tens of millions of galaxies and quasars, con¬ 
structing a 3-dimensional map spanning the nearby uni¬ 
verse to 10 billion light years. DESI will be conducted on 
the Mayall 4-meter telescope at Kitt Peak National Ob¬ 
servatory. It is supported by the Department of Energy 
Office of Science to perform Stage IV dark energy mea¬ 
surement (see [13 for definition of Stage IV-class survey). 
The expected number density of the galaxies in terms of 
the CO- moving volume is summarized in Table HI 


where (fc, /i) denotes the fiducial coordinates for the un¬ 
derlying cosmological model, and (g, v) represents the co¬ 
ordinates in the true cosmology. 

The A-P effect for bispectrum is also modeled in a 
similar way, and the resultant shape of bispectrum de¬ 
pends on five parameters, i.e., (fci, fe 2 , fca,/^ 2 )- The 
observed bispectrum is thus related to the true one given 
in Eq. (EH) through. 


H°‘’®(fcl,fc2,fc3,Ail,Ai2) = 




true \ 2 \ ^ 


j V 


X H(gi,g2,g3,!Zi,i/2) ■ (13) 

The relations between two coordinates are give by, 

qi = a{ni)ki , (14) 

and 

,, Tjtrue 

Mi P 




(15) 


z 

rig [/i^Mpc 

Vsurvey \h GpC ] 

0.6-0.8 

1.2 X 10"^ 

5.3 

0.8-1.0 

1.1 X 10“® 

7.0 

1.0-1.2 

5.4 X 10"^ 

8.3 

1.2-1.4 

3.3 X lO”"* 

9.4 

1.4-1.6 

1.5 X lO"'* 

10.1 

1.6-1.8 

5.0 X 10“® 

10.6 


TABLE I: The expected number density of galaxies Ug and 
survey volume Vurvey at each redshift bin used in the Fisher 
matrix analysis. These specific values are taken from those 
assumed in the DESI experiment. 

In order to compute the Fisher matrix, the error co- 
variance of the power spectrum and bispectrum needs to 
be evaluated. To simplify the analysis, we will ignore the 
off-diagonal components of the covariance matrices aris¬ 
ing mainly from the nonlinear mode coupling. This would 
certainly lead to an optimistic estimation of the param¬ 
eter constraints, however, it has been shown that in the 
case of power spectrum, the non-Gaussian error contri¬ 
bution to the off-diagonal components are small enough. 
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and the diagonal components of the covariance matrix 
can be approximately described by the simple Gaussian 
contribution [d^l- This would be true as long as we con¬ 
sider the quasi-linear scales at moderately high redshift 
(say, z>l and k < 0.15/iMpc“^). The parameter es¬ 
timation study with non-Gaussian covariance further re¬ 
vealed that the size of the constraints on each parameter 
is not drastically changed if we consider the multiple pa¬ 
rameter estimation ([3,14^, see also 0)- 

The Gaussian contribution to the covariance matrix 
for power spectrum is given by, 


with the matrix M given by M = (Cpp — 
CpBCgBCBp)”^. Note that with the Woodbury formula 
M = Cpp —C'ppC'pb(C'bb -I-C'bpC'ppC'pb)~^C'bpC'pp, we 
easily verify that C C~^ = C~^ C = I. 

With the full covariance matrix given above, the Fisher 
matrix combining the power spectrum and bispectrum 
becomes 


Fap 


g tgobs g^obs 

“ ^ dxa dxjj 


(23) 




1 2 


ZUk,^^)P{k) + — 

n„ 


where Np is given by, 


Np = 




survey j 2 


2(27r) 


k^AkAfi. 


(17) 


where Xa indicates the parameters for our interest, i.e., 
Da, and so on. The quantity is the signal 

vector expressed as 


(18) 


/ P°'^%k) \ 

\B°'^‘^{ki,k2,k3) ) 


(24) 


Here the survey volume is derived by calculating the co¬ 
moving shell in each redshift bin, and multiply by frac¬ 
tional factor from DESI /sky = 14,000 deg^. The Gaus¬ 
sian contribution to the covariance matrix for bispectrum 
is expressed as, 


GBcB^ = Safi SB 




Nf 


Z^,{ki,fii)P{ki) 


(19) 


Zf{k2,H2)P{k2) H- 


ZUk3,i^3)P{k3) + — 


where Nb is given by, 


Nb = 2TTkik2k3{Ak)^ (Afj,)^ 




1 2 


(27r)^ 


( 20 ) 


The pre-factor sb is set to 6 for equilateral, 2 for isosce¬ 
les, and 1 for general triangular configurations. Finally, 
the cross-covariance matrix between power spectrum and 
bispectrum is given by 


S ( 






€3 


) [ZMk) 4 

' 2Z2 [ZfPjkP + n-i] [ZfP{k 2 ) + n-^] 
Zi{ki,^i)Zi{k2,H2) 


.-1] 


( 21 ) 


cyclic 


Notice that while the covariance matrices, Cpp and 
Cbb, can be expressed as a diagonal form, the full matrix 
of C combining both power spectrum and bispectrum 
is no longer diagonal even in the Gaussian case. This 
generally requires a complicated matrix algebra, however, 
the covariance matrix C in our case can be expressed in 
a block-diagonal form as. 

Then the inverse matrix C~^ is expressed as, 

M -MCvbC^I \ 

V -CbbCbpM Cb^ + CbbCbpMCpbCb,^ ) 


with P°'°^ and B°'°^ being given by Eqs. (fl^ and (fl^ . 
Note that in the absence of the cross covariance Cpp, 
Eq. (1^ is reduced to a simplified form: 


FaP -> 

k 


^pobs(^) dP°'=%k) 

dXa dxp 


(25) 


+ E 


dB°'=^{ki,k2,k3) dB°^^{ki,k2,k3) 

6xa dxp 


Below, we will present the results in four different 
cases; 1) power spectrum only, 2) bispectrum only, 3) 
power spectrum and bispectrum, but cross covariance ig¬ 
nored [Eq. (l2^ ]. and 4) full information taking account 
of the non-vanishing cross covariance [Eq. (l23l) ]. 


III. RESULTS 


In this section, we present the results of the Eisher ma¬ 
trix forecast described in the previous section. In what 
follows, the cosmological parameters needed to compute 
the power spectrum and bispectrum are set to the Planck 
2013 concordance model [3; = 0.32, Hf, = 0.049, 

h = 0.67, As = 2.15 X 10“® and ns = 0.96. In our Fisher 
matrix, the number of free parameters to be estimated 
is six, i.e., geometric distances Da and 77“^, growth of 
structure /, galaxy biasing b and 62 , and FoG damping 
CTp. While the fiducial values of the first three parame¬ 
ters are specified by the cosmological parameters of the 
underlying cosmological model, the last three parameters 
are nuisance parameters, for which we assume that the 
observed galaxy distribution faithfully traces the mass 
distribution. That is, the fiducial values for the galaxy 
bias h and 62 are respectively set to 1 and 0 , and the 
fiducial value of Op is just given by the linear theory pre¬ 
diction. 
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FIG. 1: Signal-to-noise of power spectrum (dotted) and bis- 
pecrum (dashed) as function of redshift. Assuming a DESI- 
like experiment, the signal-to-noise ratios defined at Eqs. (I26II 
and (Hzl) are estimated at each redshift bin. The combined 
result of the power spectrum and bispectrum, taking account 
of the cross covariance, is plotted as solid curve [see Eq. (I28II 
for definition], and, taking account of the diagonals only, is 
plotted as dot-dashed curve. 


A. Signal-to-noise ratio 

Before discussing the constraints on each parameter, 
we first look at the signal-to-noise ratio S/N, and check 
the potential power of the bispectrum relative to that 
of the power spectrum. We compute the S/N for both 
the power spectrum and bispectrum at each redshift bin. 
The signal-to-noise ratio for the power spectrum and bis¬ 
pectrum, {S/N)p and {S/N)b, are respectively defined 
as. 



Similarly, we define the signal-to-noise ratio for the com¬ 
bined case: 

(I)' = 1 : 1 : ‘SC-'S. (28) 

k kiMM 

Here, the vector S is similarly defined as Eq. (1^ . but 
the quantities and are replaced with P and B 
[Eos. 151 and [TT]. 


In Eig. [U we plot three different {S/N) defined above 
and one {S/N) without the full covariance combina¬ 
tion, adopting the specific survey design of DESI. Note 
that in all cases, the cut-off wavenumber k is set to be 
k = 0.1/iMpc“^. The dotted and dash curves repre¬ 
sent the S/N for power spectrum and bispectrum, re¬ 
spectively. The signal-to-noise for bispectrum is basi¬ 
cally smaller than that for power spectrum. Neverthe¬ 
less, the combination of the bispectrum with the power 
spectrum helps to improve the signal-to-noise ratio, as 
shown by the solid curve, and the improvement becomes 
significant at lower redshift. This is because the num¬ 
ber density of galaxies becomes larger at lower redshift 
in our setup (Table |T|. Then, the shot noise contribu¬ 
tion can be suppressed at relatively larger wavenumber. 
The combined {S/N) with diagonals only is presented 
as a dot-dashed curve. Given the fact that the number 
of available configurations or modes in the bispectrum 
more rapidly increases with the wavenumber than that 
in power spectrum, {S/N)b eventually exceeds {S/N)p 
at the first redshift bin, although the cosmic variance er¬ 
ror now comes to play an important role and the total 
signal-to-noise is slightly reduced compared to the one at 
the second redshift bin. 


B. Constraints on the geometric distances 

As we mentioned in Sec. HI a precision measurement of 
the BAO scale is the key to determining the geometric 
distances, Da{z) and H~^{z), through the A-P effect. 
Although the acoustic structure imprinted on the power 
spectrum and bispectrum actually depends on cosmol¬ 
ogy, the counterpart of the acoustic oscillations is pre¬ 
cisely observed in the CMB anisotropies, and with the 
cosmological results by WMAP and Planck as prior in¬ 
formation, the BAO can be used as standard ruler. At 
/c < 0.1 hMpc~^, the galaxy power spectrum and bispec¬ 
trum are supposed to be described by the leading-order 
perturbation, where the acoustic signature is clearly vis¬ 
ible. In the presence of galaxy bias, the only uncertainty 
is the clustering amplitude, however, this does not seri¬ 
ously affect the measurement of the characteristic scales 
of BAO. 

In Fig. [21 the forecast result for the statistical errors 
of the angular diameter distance Da is shown. The frac¬ 
tional errors are presented in the lower panel. With the 
DESI-like survey, constraints on Da from the power spec¬ 
trum (dotted) and bispectrum (dashed) become compa¬ 
rable at lower redshift, 2 : < 1.2, and precision is typi¬ 
cally 1 — 2% level, although the constraining power of 
the bispectrum becomes rapidly worsen as increasing the 
redshifts. Note that if all the other parameters such as 
growth rate / and galaxy bias are known a priori, Da can 
be measured with sub-percentage accuracy. However, 
the growth rate and the linear bias are mostly degenerate 
with Da, particularly in the limit of /r —>■ 0. Combining 
both the power spectrum and bispectrum, depicted as 
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FIG. 2: Upper. The expected Icr error of the angular diame¬ 
ter distance Da{z) for DESTlike experiment. The size of the 
errorbars corresponds to the marginalized Icr error estimated 
from the Fisher matrix, combining both the power spectrum 
and bispectrum. Dotted and dash curves represent Da{z) for 
different dark energy models with equations of state parame¬ 
ter w = —1.05 and —0.95 respectively. Lower. The fractional 
errors of Da using DESI are presented. Meaning of the line 
types are the same as in Fig. [T] 

dot-dashed line, the constraint on Da will be improved 
by a factor of 1.5, thus mostly achieving a sub-percent 
precision. In the upper panel, we plot several curves for 
Da varying the dark energy equation of state parame¬ 
ter: w = —1.05 and —0.95 represented by dotted and 
dash curves respectively. The dark energy model will be 
severely constrained with the precision measurements of 
Da- 

Note that a comparable result of the constrained val¬ 
ues from the power spectrum and bispectrum is under¬ 
stood from the signal-to-noise in Fig. [TJ The signal- 
to-noise for bispectrum is higher than that of the power 
spectrum at low redshift bins, at which the number of 
observed galaxy is higher than the critical density of 
nga\ = 10“^ (/i^Mpc“^) (see Table U). Accordingly, the 
cosmological constraints on distances from bispectrum 
become better than those from power spectrum in the 
first two redshift bins, as shown in Fig.[2j Even if rigai de¬ 
creases but still exceeds 10“"^ (/i^Mpc”^), the constraints 
from bispectrum are comparable to those from power 
spectrum. Thus, combining both power spectrum and 
bispectrum provides an opportunity to improve the con¬ 
straints on Da- 

On the other hand, the determination of radial dis¬ 
tance H~^ is more difficult than that of Da- The redshift 



FIG. 3: Upper. Expected Icr errors on the comoving ra¬ 

dial distance, given by the inverse Hubble parameter 
assuming DESI-like survey. The results combining both the 
power spectrum and bispectrum are shown. Dotted and dash 
curves represent H~^{z) for different dark energy models with 
equations of state parameter w — —1.05 and —0.95 respec¬ 
tively. Lower. The fractional errors of using DESI are 
presented. Meaning of the line types are the same as in Fig. [T] 


dependence of the error on is presented in Fig. 
The dotted and dashed curves represent the constraints 
using the power spectrum and bispectrum, respectively. 
The information on H~^ is most imprinted at the limit of 
/i —> 1, at which the observed spectrum is influenced by 
RSD, mainly due to the non-linear smearing effect. Be¬ 
cause of this, the resultant statistical precision is poorer 
than that of Da by a factor 2. Comparing between the 
solid and the dot-dashed curves, the off-diagonal ele¬ 
ments of the covariance matrix do not significantly re¬ 
duce the constraining power. For refrence, in the upper 
panel, we also plot H~^ for various dark energy models 
with w = —1.05 and —0.95 represented by dotted and 
dash curves respectively. 


C. Constraints on the growth of structure 

At the scales close to the linear regime, the improved 
theoretical model of RSD successfully describes the ob¬ 
served spectrum. Still, however, nonlinear cross talk of 
the small-scale physics with large-scale anisotropic clus¬ 
tering (i.e., FoG effect) is significant, and a proper treat¬ 
ment of this cross talk needs to be incorporated. Here, we 
adopt a simple Gaussian form to phenomenologically de- 




































FIG. 4: Two-dimensional contour of the expected Icr con¬ 
straint on the parameters / and Up. The results at the redshift 
bin of 0.6 < 2 ; < 0.8 is particularly shown. The dotted con¬ 
tour represents the case of power spectrum, the dash contour 
represents the case of bispectrum, the dot-dashed contour rep¬ 
resents the case of combination with diagonals only, and the 
solid curve represents the case of the full covariance combina¬ 
tion. 

scribe the FoG effect [see Eqs. O and mi This would 
certainly remedy the flaw in the perturbative modeling 
of the RSD, however the strength of the damping char¬ 
acterized by fjp is now a free parameter, which needs to 
be determined by observations. 

Fig. |T] shows the impact of this uncertainty in con¬ 
straining the growth of structure. Here, the plotted result 
is the two-dimensional contour of the expected Icr con¬ 
straint, which is estimated at the specific redshift bin, 
0.6 < z < 0.8. The parameter tTp is significantly corre¬ 
lated with the linear growth rate /, and the analysis using 
the power spectrum alone, depicted as dotted contour, 
exhibits a strong parameter degeneracy. This is partly 
due to the small cutoff wavenumber, fc^ax = 0.1 /iMpc“^, 
below which the damping term, expressed as function of 
, is monotonically varied along the direction fi, and 
the behavior looks very similar when we vary /. The 
situation almost remains unchanged even if we use the 
bispectrum, and the strong degeneracy between / and 
CTp is observed (dashed contour). Interestingly, however, 
the direction of the parameter degeneracy differs from 
that of the power spectrum case. Accordingly, the com¬ 
bination of the power spectrum and bispectrum improves 
the constraint on /. The dot-dashed contour represents 
the case combination with diagonals only. 

Figs.[S]and|n]respectively show the marginalized results 


of the fractional errors on / and Up, plotted as a function 
of redshift. Combining power spectrum and bispectrum, 
the constraint on / becomes tighter at all redshift bins. 
In particular, at lower redshift bins, statistical power of 
the bispectrum is enhanced, and the relative impact on 
the constraint on the parameter cTp eventually becomes 
stronger (see bottom panel of Fig. [SI). As a result, com¬ 
pared to the power spectrum results, the combined result 
of the constraint on / is improved by a factor of 2. Note 
that similar to the constraints on the geometric distance, 
the effect of the off-diagonal components of the error co- 
variance are insignificant in the case of / and Up, and it 
does not appreciably change the results. In the upper 
panel, we plot / for different theoretical models. The 
black dotted and dash curves represent / for dark energy 
models with w = —1.05 and —0.95 respectively. It is in¬ 
teresting to note that we are able to constrain the dark 
energy from the measured / alone. Also, the red long 
dash curve represents / of DGP model [48 - l^ . demon¬ 
strating the outperformance of the stage-IV class survey. 

Finally, to see the impact of the FoG effect on the 
measurement of the growth rate, we plot in Fig. |7| the 
marginalized constraints of / in the case when we know 
the FoG effect a priori. As it is expected, the parameter 
/ is better constrained, and the precision reaches at ~ 
2% level. One notable and interesting point is that the 
power spectrum always gives a tighter constraint, and 
combination of power spectrum and bispectrum does not 
improve much the constraint. Thus, the benefit of using 
the bispectrum is substantially reduced. 


D. Constraints on the galaxy bias parameters 

The combination of power spectrum and bispectrum 
has been frequently used in the literature to constrain the 
bias parameters (e.g., (s^, [s^ for recent results). With 
the survey like DESI, the constraint on these parameter 
will become much more tigher. This is true even if we si¬ 
multaneously estimate the expansion history and growth 
of structure. 

Fig. [ 8 ] shows the forecast result of the statistical errors 
on galaxy bias parameters. The upper panel presents the 
constraint on the linear bias parameter b. The linear bias 
is tightly constrained by the power spectrum depicted as 
dotted line, and no benefit to use the bispectrum is found 
(dashed), although the combination of power spectrum 
and bispectrum still improves the constraint by a factor 
of 2 (solid). This is because of another bias parameter 
62 , which coherently boosts the amplitude of bispectrum. 
Also, one important assumption here is that the cosmo¬ 
logical parameters are known a priori from the GMB ob¬ 
servations, including the normalization of linear power 
spectrum. In other words, if the specific cosmological 
models are not given, the linear bias information is not 
extracted from the power spectrum. 

On the other hand, the bottom panel of Fig. [5] shows 
the constraint on the nonlinear bias parameter 62 - The 




9 



FIG. 5: Upper. Expected la constraint on the linear growth 
rate f{z), assuming DESI. The results are estimated using 
both the power spectrum and bispectrum. Dotted and dash 
curves represent f{z) with w = —1.05 and —0.95 respec¬ 
tively, and red long dash curve represents f{z) with DGP 
model.iouier; The fractional errors of /. Meaning of the line 
types are the same as in Fig. [T] 


power spectrum alone cannot constrain 62 , because the 
parameter 62 does not appear in the expression of power 
spectrum, relevant at the large scales of our interest (see 
Sec. Ill All . Even using the bispectrum, 62 is poorly con¬ 
strained. This is due to the degeneracy between the bias 
parameters b and 62 - In this respect, the combination of 
the power spectrum and bispectrum is quite essential to 
break the degeneracy. With the help of the power spec¬ 
trum information, b is measured precisely, and this will 
lead to a substantial improvement of the constraint on 
62 , as shown by the solid line. 


IV. SUMMARY 

The bispectrum has been recognized as a powerful tool 
to probe non-Gaussian initial conditions, and to enhance 
the constraints on BAO distance measurements. In this 
paper, we highlighted the benefit to improve the simulta¬ 
neous constraint on the growth of structure, and geomet¬ 
ric distances through the RSD and A-P effect. The power 
spectrum in redshift space suffers from nonlinear cross 
talk between small- and large-scale clustering, and the 
clustering amplitude is signihcantly suppressed at large 
scales. While the impact of this FoG effect is less signifi¬ 
cant at the scales close to the linear regime, it is difficult 



FIG. 6 : Upper. The expected error on nuisance parameter 
characterizing the damping scales of the FoG effect, ap{z). 
The resnlts are presented assnming DESI-like experiment. 
Lower. The fractional errors of ap{z). Meaning of the line 
types are the same as in Fig. [T] 



z 


FIG. 7; Fractional errors of the linear growth rate / in the 
case when we know a priori the FoG effect (i.e., ap). Meaning 
of the line types is the same as in Fig. [T] 
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FIG. 8 : Expected error on galaxy bias parameters. Up¬ 
per panel shows the fractional errors on linear bias 6 , while 
the lower panel presents the errors on second-order bias 62 . 
Meaning of the line types is the same as in Fig. [2] 

to break degeneracy between the distortions induced by 
coherent motions and systematics arising from the FoG 
effect. Although this would be certainly resolved with a 
better understanding of the nonlinear cross talk based on 


a more elaborate theoretical model, we proposed an al¬ 
ternative method to improve the cosmological constraints 
by using both the power spectrum and bispectrum. As¬ 
suming a DESI-like experiment, Fisher matrix analysis 
suggests that the analysis using the power spectrum data 
alone shows a strong parameter degeneracy between the 
growth of structure (/) and ap as a nuisance parameter 
characterizing the damping scale of the FoG effect. In 
hat case the estimated error on / is about 5% at the 
most sensitive redshift bin. We found that this parame¬ 
ter degeneracy can be broken when combining both the 
power spectrum and bispectrum. As a result, the con¬ 
straint will be improved by a factor of two, which satisfies 
the primary science goal to distinguish the possibilities 
to explain the cosmic acceleration, i.e., the dark energy 
and modification of gravity on cosmological scales. 

Note finally that these results certainly depend on the 
model of the FoG effect. This is especially the case for 
the bispectrum. While we adopted a simple Gaussian 
ansatz to describe the FoG damping, a proper way to 
characterize the FoG effect in bispectrum is not yet fully 
understood. Different assumption or prescription of the 
FoG effect is possible, and may lead to a quantitatively 
different result. We leave this issue for further work. 
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